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Abstract. A symmetric nx n matrix A is called completely positive (CP) if 
there exists an entrywise nonnegative matrix B £ M. nxm that A = BB T . 
Given a nonnegative P-matrix whose entries are partially given, the CP- 
completion problem is to decide whether we can assign values to the ungiven 
entries such that the completed matrix is CP. The main results of this pa- 
per are as follows: 1) We first present an SDP algorithm to solve the CP- 
complction problem. 2) When all the diagonal entries of the P-matrix are 
given, the algorithm can give a certificate for the non-CP-completcncss, and 
is almost always guaranteed to get a CP-completion if it is CP-complctable. 
3) When the diagonal entries of the P-matrix are partially given, if the corre- 
sponding maximal principal submatrix is not CP-complctable, then a certifi- 
cate can be obtained; if the maximal principal submatrix is CP-complctable, 
then the algorithm is almost always guaranteed to get a CP-complction. In 
particular, if the P-matrix has only one diagonal given, then it can always be 
completed to a CP-matrix. In all our computational experiments, we always 
get a CP-complction if it exists. 



1. Introduction 

A real nx n symmetric matrix A is called completely positive (CP) if there exist 
nonnegative vectors u\, ■ ■ ■ , u m such that 

(1.1) A = uiuj H hfi ra !i^, 

where the number m is called the length of the factorization (1.1). The smallest pos- 
sible m is called the CP-rank of A. If A is CP, we call (1.1) the C P -decomposition 
of A. Equivalently, A is CP if it could be factorized as A = BB T , where B e R nxm 
is an entrywise nonnegative matrix. Hence a CP-matrix is not only positive semi- 
definite but also nonnegative. 

CP-matrices have wide applications in mixed binary quadratic programming [7], 
approximating stability numbers [10], the MaxClique problem [29,33], single qua- 
dratic constraint problem [30], standard quadratic optimization [3], general qua- 
dratic programming [31], and so on. Some NP-hard problems can also be refor- 
mulated as linear optimization problems over the cone of CP-matrices [16,21,23]. 
These important applications motivate researchers to study whether a matrix is 
CP or not. The interested reader may refer to [1,2,4-6,8, 11-13] for more details. 

An interesting problem related to the CP-decomposition is the CP-completion 
problem. Let E = {(h,ji),--- , | 1 < ik < 3k < n, k = 1, ■ ■ ■ ,1}. A real 

n x n symmetric partial matrix is called an E-matrix, if its entries are given for 
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all € E, while all the other entries are ungiven. The CP-completion problem 
is to study whether we can assign values to the ungiven entries of the P-matrix 
such that the completed matrix is CP. As a symmetric matrix can be identified 
by a vector that consists of all its upper triangular entries, we take the given upper 
triangular entries to form the identifying vector of an P-matrix as follows: 

a G R E , 

where M. E denotes the space of vectors that are indexed by (i, j) € E. It is obvious 
that an P-matrix with an identifying vector a S M. E is CP-completable if and only if 
there exists a CP-matrix A such that a = A\e, where A\e denotes the P-truncated 
vector of the matrix A. We say an P-matrix is non-CP -completable if it could not 
be completed to a CP-matrix. 

For example, consider the P-matrix [1, Example 2.23]: 

2 3 * 

3 6 3 
3 6 3 
* 3 2 



(1.2) 



where * denotes the ungiven entries of the P-matrix throughout this paper. The 
index set of (1.2) is P = {(1, 1), (1, 2), (1, 3), (2, 2), (2, 3), (2, 4), (3, 3), (3, 4), (4, 4)} 
and the corresponding identifying vector is a = (2, 3, 0, 6, 3, 0, 6, 3, 2) T . Note that 
for any nonnegative * in the (l,4)-th position, the determinant of (1.2) is det = 
-27(* + 1) < 0, so it is not CP-completable. 
If we consider the P-matrix: 

"ill 

(1.3) 1 1 1 

1 1 * 

Then the index set is P = {(1, 1), (1, 2), (1, 3), (2, 2), (2, 3)} and the corresponding 
identifying vector isa=(l,l,l,l,l). Since 

T 



A = 




is completely positive and a = A\e, (1-3) is CP-completable. 

Note that if the diagonal of a CP-matrix is zero, then all the off-diagonal entries 
on this row and column will be zeros too. Without loss of generality, we assume 
that the P-matrix we wish to analyze is nonnegative and all its given diagonals arc 
strictly positive. 

There are some results on CP-completion, which are related to the underlying 
graph of the P-matrix with the special structure. An P-matrix is called a partial 
CP-matrix if all its fully specified principal submatrices are CP. Every partial CP- 
matrix realization of a connected graph G is CP-completable if and only if G is a 
block-clique graph [1,14]. If G is the specification graph of the P-matrix, we call the 
P-matrix a G -partial matrix. For a partial GP-matrix whose specification graph 
contains cycles, under some additional cycle conditions, the G-partial CP-matrix 
is CP-complctablc if and only if all the blocks of G are cliques or cycles [15]. 

The above CP-completion results are focused on the special P-matrix that have 
all diagonals given or the corresponding specification graph satisfies some struc- 
tures. In this paper, we will study the CP-completion problem for the general 
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P-matrix. We will discuss the problem in three cases regarding whether the diago- 
nals of the i?-matrix are given or not. The first case is that all diagonals are given, 
and the other two cases are that all diagonals are ungiven and the diagonals are par- 
tially given. We will explore if an ©-matrix can not be completed to a CP-matrix, 
how can we get a certificate for that? If an ©-matrix is CP-complctable, how can 
we get a CP-completion and give a CP-decomposition for this completion? 

The paper is organized as follows. In Section 2, we analyze the CP-completion 
problem for the P-matrix with all diagonals given. We show that the CP-completion 
problem can be formulated as an £ -truncated A-moment problem (£-TAMP) with 
some special £ and A. We present an SDP algorithm for the £-TAMP; some prop- 
erties of the algorithm are obtained. In Section 3, we discuss the CP-completion 
problem for the P-matrix with the diagonals ungiven and partially given, respec- 
tively. Some experimental results are given in Section 4. Finally, we conclude the 
paper in Section 5. 

2. AN SDP ALGORITHM FOR CP-COMPLETION 

Recently, Nie [24] presented a new efficient SDP algorithm for the £ -truncated A- 
moment problem (f -TAMP), which is a generalization of the scmidcfinitc approach 
proposed in [18] for the truncated A-moment problem (TAMP). In this section, we 
will show how to formulate the CP-completion problem as an £-TAMP with a 
special £ and a compact set A. We then present an SDP algorithm for it, which is 
a special case of Nie's algorithm. Some basic properties of the algorithm will also 
be studied. 

2.1. Formulation as f-TAMPs. We begin with a sufficient and necessary con- 
dition for an P-matrix to be CP-completable. Let 

(2.1) A = {x G R" : h{x) = 0, g(x) > 0} 

be a standard simplex, where h = Xi + ••• + x n — 1 and g = (xi, • • • , x n ) are 
two tuples of polynomials. It is easy to see that an P-matrix with the identifying 
vector a G M. E is CP-completable if and only if there exist nonnegative vectors 
vi , • • • , v m G A and p\ , • • • , p m > such that 

(2.2) A := pivivf H h p m v m v^ and a.= A\ E - 

Denote the set of nonnegative integers by N. For any a := (ai, . . . , a„) T G N n , we 
define \a\ := a± H + a n . Let 

(2.3) £ = {a G N" : \a\ = 2,<p{a) G P}, 

where ip(ot) is the following one-to-one mapping from the index vector set {a G 
N n : \a\ = 2} to the point pair set ■ 1 < i < j < n}: 

a= (0,--- ,0,1 4 ,0,--- ,0,1„0--- ,0) T ^(i,j), 

or 

a = (0,-- - ,0,24,0,- •• ,0) T ^ 

with li and % meaning the i-th entry of the vector is 1 and 2 respectively. We say 
£ is the index set of P. For example, if n = 3 and P = {(1, 2), (2, 2), (2, 3)}, then 
£ = {(1, 1, 0), (0, 2, 0), (0, 1, 1)}. The degree of £ is defined as deg(£) := max{\a\ : 
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a £ £ }. Thus deg(£) = 2 for any E. Hence, with a proper term order, we could 
rewrite the identifying vector a as 



where a Q denotes the entry indexed by a and M. £ denotes the space of real vectors 
indexed by elements in £. We call a an £ -truncated moment sequence (f-tms). 

Let A be as in (2.1). The £-truncated A-moment problem (£-TAMP) studies 
whether or not a given £-tms a admits a A-measure p, i.e., a nonnegative Borel 
measure \x supported in A such that a a = J A x a dp for any a £ £, where x a := 
x" 1 ■ ■ ■ a;"" . A measure p, satisfying the above condition is called a A-representing 
measure for a. Furthermore, we say a measure is finitely atomic if its support is 
a finite set, and is m-atomic if its support consists of at most m distinct points. 
In other words, if an £-tms a admits a m-atomic A-measure, then there exist 
Vi , . . . , v m £ A and p\ , . . . , p m > such that a = p\ [i>i]g + • • • + p r [v m )s , where 
[ v i]s '■— { v f)a£S f° r i — 1, • ■ • , m. Hence by (2.2), an ^-matrix with the identifying 
vector a 6 M £ is CP-complctable if and only if a admits a m-atomic A-measure, 
i.e., 

(2.5) a = p\[vi]e H h Pm[v m }s, 

where Vi £ A and pi > (i — l,...,m). Now we have formulated the CP- 
completion problem as an £-TAMP with £ and A given by (2.3) and (2.1) respec- 
tively. 

2.2. A semidefinite algorithm. We first review some basic conceptions. Denote 
M.[x]e := span{x a : a £ £}. R[x]s is said to be A- full if there exists a polynomial 
p(x) £ R[x]s such that p{x)\a > 0, which implies that p(x) is always positive on 
A [28]. Denote N% := {a € N n : \a\ < d} and R[x] d := span{x a : a G N%}. 
The dimension of R[x] d is s(d) = C% +d = [22]. For a tms z £ and 

q £ M\x\2ki a symmetric matrix Lq k \z) is called the k-th localizing matrix of q 
generated by z if it is linear in z and satisfies the Riesz function. The order of 
Lq k \z) is s(k — \deg(q)/2~\), where \t\ is the smallest integer that is not smaller 
than ( £ M [24]. Especially, L±(z) is called the k-th order moment matrix and 
denoted by Mk(z). 

Let go — 1 and g n +i := 1 — IMH- Since the set A defined by (2.1) satisfies 
A C B(0, 1) := {x £ 1™ : < 1}, it can be equivalcntly described as A = {x £ 
R n : h = 0, gs '■= (9o,9x, ■ ■ ■ 7 5n,5n+i) > 0}. If n = 2 and k = 2, then the k-th 
localizing matrices of polynomials h and gs could be computed as follows: 



z (0,0) z (l,0) 2 (0,1) z (2.0) Z (1S) z (0,2) 

z (lfi) z {2,0) Z {1,1) z (3,0) Z (2S) z (l,2) 

M2 f z \ = Z (0-l) z (l,l) z (0,2) z {2.\) 2(1,2) Z(0,3) 

2(2,0) 2(3,0) 2(2,1) 2(4,0) 2(3,1) 2( 2) 2) 

2(1,1) 2(2,1) 2(i,2) 2(3,1) Z( 2 ,2) Z(l,3) 

2( ,2) 2(i,2) 2( ,3) Z( 2 ,2) 2(i, 3 ) 2( , 4 ) 

2(1,0) + 2(0,1) - 2(0,0) 2(2,0) + 2(1,1) - 2(i,o) 2(1,1) + Z( ,2) ~ ^(0,1) 



(2.4) 



a = (a a ) ae£ £ R £ , 




2(2,0) + 2(1,1) - 2(i,o) 2(3,o) + 2(2,1) ~ 2(2,0) 2( 2 ,1) + 2( 12 ) 
2(1,1) + 2(0,2) — 2(o,l) 2(2,1) + 2(1,2) — Z (l,l) Z (l, 2 ) + z (0,3) 



2(0,2) 
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2(1,0) 
2(2,0) 
2(1,1) 



z (2,0) 2(1,1) 
2(3,0) 2(2,1) 
2(2,1) «(1,2) 




2(0,1) 2( M ) Z( 02 ) 
2(1,1) 2(2,1) 2(1,2) 
2(0,2) 2(1,2) 2(0,3) 



2(0,0) _ 2(2,0) _ 2(o,2) 2( 10 ) - Z( 30 ) - Z( 12 ) Z( ,l) 
L l- x j-xl (2) = 2(i, 0) - Z( 3 ,o) - 2(i,2) 2(2,0) - 2(4,0) ~ 2(2,2) Z(l,l) 
2(o,i) - 2(2,1) - 2(o,3) 2(i,i) ~ 2(3,1) - 2(1,3) 2(o,2) 



2(2,1) - 2(0,3) 
2(3,1) - 2(1,3) 
2(2,2) - 2(0,4) 



It is easy to know that if z G ffi 2fc admits a A-mcasure //, then 



(2.6) 



L^(z) = 0,L^(z)hO, (7=0,1,..., 



n + 1), 



where L 9 . (z) >r denotes that Lg.'(z) is positive semidefinite. Hence (2.6) is 
necessary for z to admit a A-measure. However, the converse is not always true. 
We say z is flat if it satisfies not only (2.6) but also the following rank condition 



Curto and Fialkow proved in [9] that if a tms z 6 R N 2fc is flat, then z admits a 
unique A-measure, which is rankAffe(z)-atomic. 

Given two tms' y G IR Nd and z G R N = , we say z is an extension of y, if d < e 
and y a = z a for all a G N^- In general, we denote by z\g the subvector of z, whose 
entries are indexed by a £ £. For convenience, we denote by z^ the subvector 
If z is flat and extends y, we say z is a /?a£ extension of y. Indeed, it is 
shown in [24, Proposition 3.3] that an £-tms a G R £ admits a A-measure if and 
only if it could be extendable to a flat tms z G R N 2<= for some k. Hence based on 
(2.5), determining whether an E- matrix with an identifying vector a G M. £ is CP- 
completable could be transformed to investigating whether a has a flat extension. 

We now start to solve the CP-completion problem which has been formulated 
as the £-TAMP. Let d > deg(8) — 2 and R £ M.[x]d- Denote the coefficient vector 
of the polynomial R in the graded lexicographical ordering by (R a )aeN™- Consider 
the following linear moment optimization problem 



where T^(A) = {z £ K. N d : meas(z, A) ^ 0} with meas(z, A) being the set of all 
A- measures admitted by z. 

Note that A C B(0, 1) is compact. If M.[x]e is A-full, then the feasible set of 
(2.8) is compact convex and (2.8) has a minimizcr for all R. If M[x]s is not A-full, 
it needs to choose R which is positive definite on A to guarantee that (2.8) has a 
minimizcr. Therefore, to get a minimizer of (2.8), it is enough to solve (2.8) for 
a generic positive definite R, no matter whether M.[x]s is A-full or not. Thus, we 
choose a generic R £ where is the set of all sum of squares polynomials 

in n variables with degree d. 

Since T^(A) is quite difficult to describe, we relax it by (2.6) to the following set 



(2.9) T k (h,g B ) :={zeK N - | hf (z) = 0, Lf}{z) h 0,j = 0, 1, . . . , n + l} , 



(2.7) 



rankAffc_i(z) = rankM^(z). 



(2.8) 



min R 

s.t. z\ e = a, z g T d (A), 
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where k > d/2 is an integer. Then the k-th order semidefinite relaxation of (2.8) 
can be stated as 

min RaZ a 
(2.10) (SDR) k : z «6Nj fc 

s.t. z\e = a,z G T k (h,g B )- 

Now we can use the following hierarchy SDP algorithm to solve the £-TAMP. As 
mentioned above, we start with an even d > deg{£) = 2 and an order k > d/2. 

Algorithm 2.1. The SDP algorithm for the £-TAMP 
Step 0: Choose a generic R G £ n ,d> an d let k := d/2. 

Step 1: Solve (2.10). If (2.10) is infeasible, then a doesn't admit A-measures, 
and stop. Otherwise, compute a minimizer z* ,k . Let t := 1. 

Step 2: Let w := z*' k \2t- If the rank condition (2.7) is not satisfied, go to Step 

4. 

Step 3: Compute the finitely atomic measure [i admitted by w: 

H = Ai<5(wi) H h X m S(um), 

where m = rankA/ t (u>), Ui G A, Xi > 0, and 5(iii) is the Dirac measure supported 
on the point Ui (i = 1, • • • , m). Stop. 

Step 4: If t < k, set 4 := t + 1 and go to Step 2; otherwise, set /c := fc + 1 and 
go to Step 1. 

Remark 2.2. Algorithm 2.1 is a special case of Algorithm 4.2 in [24]. Denote 
[x\d ■= (x a ) a £f$ri. We chose R = [^]J / / 2 ^ T ^Nd/2 in (2-10) where J is a random 
square matrix obeying Gaussian distribution. We checked the rank condition (2.7) 
numerically with the help of the singular value decomposition [17]. The rank of a 
matrix is evaluated as the number of its singular values which are greater than or 
equal to 10~ 6 . We used the method in [19] to get a m-atomic A-measure for w. 

2.3. Properties of Algorithm 2.1. We first present some of the main results of 
Algorithm 2.1 as follows (cf. [24, Section 5]). 

Theorem 2.3. 1) If & doesn't admit A-measures andW.[x]£ is A-full, then (SDR)^ 
is infeasible for all k big enough, which gives a certificate for the nonexistence of 
representing measures. 2) If & admits a A-measure, then for almost all generated 
R: i) We could asymptotically get a flat extension of a by solving the hierarchy 
{{SDR)k}^? =1 . ii) The obtained flat extension admits a m-atomic A-measure with 
m < \£\, where \£\ denotes the cardinality of the set £. 

Remark 2.4. Under some general conditions, which is almost sufficient and neces- 
sary, we could get a flat extension of a by solving (SDR) k for some k [25-27], which 
implies that the SDP algorithm has finite convergence, and this always occur in our 
numerical experiments. 

From the above, we know that if an .E-matrix with the identifying vector a G 
M. £ could be completed to a CP-matrix, then we could asymptotically get a flat 
extension of a for almost all R G S n ,d by Algorithm 2.1. Moreover, it could likely 
be obtained within finitely many steps. After getting a flat extension of a, we could 
get a m-atomic A-measure for a. Then we obtain a CP-completion matrix for the 
P-matrix, and a CP-decomposition is also produced. 

Note that under the condition that M.[x]s is A-full, Algorithm 2.1 could give a 
certificate for the nonexistence of CP-complctcncss. However, when M.[x]s is not 
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A-full and a doesn't admit A-measures, it is not clear whether there exists a k 
such that (2.10) is infeasible. So an interesting question is under which condition 
K^Jg is A-full? In the following, we will give a sufficient and necessary condition 
for R[x]e to be A-full. 

Proposition 2.5. Given E = ■ ■■ , | 1 < %k < jk <n,k = !,■■■ ,1}. 

Let £ be the index set of E and A = {i6 R™ : x\ + • • • + x n = 1}. Then R-[x]g is 
A-full if and only if < i < n} C E. 

Proof. We first prove the sufficient condition. If {(i,i),l < i < n} C E, then 
{(2,0,- ■■ ,0),(0,2,-- - ,0),--- ,(0,0, ••• ,2)} C £, so we have x\ £ R[x} £ for all 
1 < i < n. Hence for any x £ A, there exists a polynomial p(x) = J27=i x 1 e ^Mf 
such that p(x) > 0, i.e. p(x)\a > 0. Thus R[x]s is A-full. 

Wc prove the necessary condition by contradiction. Suppose there exists some 
io £ {1, • • • , n} such that (io, io) ^ E. Then for any polynomial p(x) £ WL[x]e, p(x) 
is a linear combination of all the monomials of degree 2 except the monomial xf . 
Let c = (0, • • • , 0, lj , 0, • • • , 0) T £ A be a constant vector, then p(c) = holds for 
all p(x) £ R[x]e. Hence, there is no polynomial p(x) £ R[x]s such that p(x)\a > 0. 
This contradicts the fact that M.[x]e is A-full. The proof is completed. □ 

Remark 2.6. Proposition 2.5 shows that R[x]f is A-full if and only if all the diag- 
onals of the .E-matrix are given. So, when the diagonals are all given, Algorithm 
2.1 could determine whether an i?-matrix could be completed to a CP- matrix 
or not. Moreover, if it is CP-completable, not only a CP-completion but also a 
CP-decomposition could be obtained by Algorithm 2.1. In fact, Algorithm 2.1 
computationally works even if there are ungiven diagonals in the P-matrix. 

3. CP-COMPLETION WITH UNGIVEN DIAGONALS 

In this section, we will study the CP-completion of the P-matrix under the other 
two cases. One is that all the diagonals of the P-matrix are ungiven, and the other 
is that the diagonals are partially given. Some properties of the P-matrix under 
these two cases will also be investigated. 

3.1. All diagonals are ungiven. For a nonnegative P-matrix with all the diag- 
onals ungiven, the following proposition shows that it can always be completed to 
a CP-matrix. 

Proposition 3.1. Suppose all the diagonals of the nonnegative E-matrix are un- 
given. Then the E-matrix can always be completed to a C P -matrix. 

Proof. Without loss of generality, we suppose the P-matrix has the following form 

* A 12 ■■■ A ln 
A 2 i * ■■■ A 2n 

A n l A n2 

where Aij(\ < i < j < n) are all given and nonnegative. Let 
(3.2) A= ]T A JJ (e 1 +e,)(e 1 + e,) T , 

l<i<j<n 

where e^ denotes the i-th unit vector throughout the paper. It is easy to check that 
A is a CP-completion of (3.1). The proof is completed. □ 
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Remark 3.2. By Proposition 3.1, we know that the P-matrix with all diagonals 
ungivcn is always CP completable. Algorithm 2.1 is guaranteed almost always to 
get a CP-completion if it exists. In our numerical experimentation, we always get 
it in fintcly many steps. Moreover, it typically gives a CP-decomposition whose 
length is much smaller than the length in the proof. 

3.2. Diagonals are partially given. In this subsection, we will study the CP- 
completion problem for the case that some of the diagonals of the P-matrix are 
given, while the other diagonals arc ungivcn. The submatrix that consists of all 
the rows and columns whose diagonals are given is called the maximal principal 
submatrix of the P-matrix . We have the following result. 

Proposition 3.3. Suppose the diagonals of the E-matrix are partially given. If the 
maximal principal submatrix of the E-matrix can not be completed to a C P -matrix, 
then the E-matrix can not be completed to a CP-matrix either. 

Proof. We prove by contradiction. Assume the P-matrix with the identifying vector 
a e R £ is CP-completable. By (2.5), we know that a admits a A-measurc. Let 
the corresponding £'-tms for the maximal principal submatrix be a'. It is easy 
to know that a is an extension of a' without considering the orders of these two 
matrices. Since a admits a A-measure fi, so does a'. Then it follows from (2.5) 
that the maximal principal submatrix with the identifying vector a' is also CP- 
complctable, which leads to a contradiction. The proof is completed. □ 

Remark 3.4. The converse of Proposition 3.3 is not necessarily true. For example, 
consider the nonnegativc P-matrix 



Since the determinant of this matrix is equal to —1 for any nonnegative * in (1, 3)-th 
position, it is not positive semidefmitc. So it can not be completed to a CP-matrix. 
However, its maximal principal submatrix 



is a CP-matrix. 

Though the P-matrix (3.3) is not CP-completable, we could show that there 
always exists a sequence of CP-completable P-matrices whose identifying vectors 
converge to the identifying vector of the P-matrix. We first give the following 
definition. 

Definition 3.5. A nonnegative P-matrix is called partial CP-complctable, if its 
maximal principal submatrix is CP-complctablc. 

Note that this definition is different from the definition of the partial CP-matrix 
in Section 1. As we know, if a matrix P is positive semidefmitc, then P + el is 
positive definite for any e > 0, where / is the identity matrix. Hence we can regard 
the positive semidefmitc matrix as a limit of some sequence of positive definite 
matrices. The converse is also true, i.e. the limit of any sequence of positive definite 
(or positive semidefmitc) matrices is a positive semidefmitc matrix. Fortunately, an 
analogous statement for partial CP-complctable matrix also holds true. We show 
it in the following theorem. 



(3.3) 



1 
1 
2 



1 2 
1 3 

3 * 



1 1 
1 1 
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Theorem 3.6. Suppose the nonnegative E-matrix with the identifying vector a 6 
M. £ is partial CP-completable. Then there exists a sequence of CP-completable E- 
matrices whose identifying vectors converge to a, no matter whether the E-matrix 
is CP-completable or not. 

Proof. If the P-matrix is CP-completable, then there exists a CP-matrix A such 
that a = A\e- Choosing A^ = A for all k = 1, 2, • • • , we know the theorem holds 
true. 

For the case that the P-matrix is not CP-completable, without loss of generality, 
we assume it has the following form 

A hn 



(3.4) 



A 1 : 

An—l,n 
A n i ■■■ Ann-i * 



where the maximal principal submatrix A' is CP-completablc and the off-diagonal 
entries in n-th row and column are all given and nonnegative. Consider the following 
sequence of P-matrices: 

Ai, n 



(3.5) 



A' + e k I n -i ■ 

A n -1, 

A n i ■ ■ ■ Ann-i * 



fc = l,2,--- 



where I n —i is the identity matrix of order n—1 and < Sk — > as k — > +oo. Hence 
the sequence of identifying vectors of (3.5) converges to a as k — > +oo. Let A' be 
a CP-completion matrix of A' , and 

A 7 



A, 



1 



. . A . A. 



nxn \<i< n -l 



where is the zero vector and b is any given nonnegative constant. Then Ak is 
a CP-completion of (3.5), which implies (3.5) is CP-completable. The proof is 
completed. □ 

Remark 3.7. By Theorem 3.6, we know that if an P-matrix is partial CP-complctablc 
but not CP-complctablc, then it must be on the boundary of the set of CP- 
complctable P-matrices, which implies that the set of CP-completable P-matriccs 
is not closed. 

Especially, when the P-matrix has only one diagonal given, we have the following 
result. 

Proposition 3.8. Suppose only one diagonal of the nonnegative E-matrix is given 
and positive. Then the E-matrix can always be completed to a CP-matrix. 

Proof. Without loss of generality, we suppose n > 2 and the P-matrix has the 
following form 

" An A 12 ■■■ A ln ~ 
A21 * ■ ■ ■ A 2n 



(3.6) 



A n l A n 2 
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where An > and all the off-diagonal entries Ay (i ^ j) are given and nonnegative. 
Let 



.4 



E ( 

l=i<j<n 

+ E 

2<i<j<n 



An 



n - 1 



1 



A ijGj f 



Aijiei + ej)(e l + ej 



Then, A is a CP-completion of (3.6). The proof is completed. 



□ 



Remark 3.9. Note that: 1) If the P-matrix is CP-completable and its diagonals 
are partially given, we could still use Algorithm 2.1 to give a CP-decomposition 
whose length is relatively small. 2) If the P-matrix is not partial CP-complctablc, 
a certificate for the nonexistence (the (2.10) is infeasible) of CP-completion for the 
maximal principle submatrix could be obtained by applying Algorithm 2.1, which 
also implies the non- CP-completeness of the P-matrix by Proposition 3.3. 3) If 
the P-matrix is partial CP-completable, then Algorithm 2.1 is guaranteed almost 
always to get a CP-completion if it exists. In our numerical experimentation, we 
always get it in finitely many steps. 



4. Numerical experiments 

In this section, we carried out some numerical experiments by using software 
GloptiPoly 3 [20] and SeDuMi [32]. We chose d = 4 and k = 2 in Step of our 
Algorithm 2.1. 

Example 4.1. Consider the nonnegative P-matrix with all the diagonals given 
(cf. [1, Exercise 2.57]): 



(4.1) 



where b is a given nonnegative constant. It is easy to check that (4.1) is CP- 
completable if and only if b > 3. Without loss of generality, we take b = 3. Clearly, 
P = {(1, 1), (1,2), (1,3), (2,2), (2,3), (2,4), (3,3), (3,4), (4,4)} and is 



2 




1 




1 





































1 









2 




1 




1 




















1 







1 









1 









2 




1 
































1 









1 




2 



We applied Algorithm 2.1 to the corresponding £-tms a = (3, 3, 0, 6, 3, 0, 6, 3, 3) T , 
and obtained the following CP-completion 



.4 



3 
3 


0.0000 



0.0000 

3 
3 



3 

E Pi u i u T> 



where the support points m and the corresponding A-rcprcsenting measure weights 
Pi are given in Table 1. Thus (4.1) with b — 3 is CP-completable, and a CP- 
decomposition is also obtained as above. 
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i 




Pi 


1 


(0.5000, 0.5000, 0.0000, 0.0000) T 


12.0000 


2 


(0.0000, 0.0000, 0.5000, 0.5000) T 


12.0000 


3 


(0.0000, 0.5000, 0.5000, 0.0000) T 


12.0000 



Table 1. The points itj and weights pi for A in Example 4.1. 



Example 4.2. Consider the nonncgative P-matrix with all the diagonals ungivcn: 



(4.2) 



* 4 

4 * 
1 
2 
2 



1 


* 

1 1 

3 2 



It follows from Proposition 3.1 that (4.2) can be completed to a CP-matrix. The 
completion matrix could be computed as (3.2). We also applied Algorithm 2.1 and 
obtained the following CP-completion matrix 



A = 



3.5379 
4 
1 
2 
2 



4 

8.4724 

1 

3 



1 



2.0113 
1 
2 



2 
1 
1 

1.5607 
1 



2 
3 
2 
1 

3.2138 



E 

i=l 



p i u 1 u l 



where itj and are shown in Table 2. 



i 


u» 


Pi 


1 


(0.3318, 0.0000, 0.2080, 0.3222, 0.1380) T 


11.6931 


2 


(0.3077, 0.0000, 0.3198, 0.1908, 0.1819) T 


1.9632 


3 


(0.0000, 0.0000, 0.4420, 0.0327, 0.5253) T 


6.6756 


4 


(0.3711, 0.3459, 0.0000, 0.2329, 0.0499) T 


3.8529 


5 


(0.2315, 0.5291, 0.0000, 0.0455, 0.1937) T 


28.6107 



Table 2. The points itj and weights pi for A in Example 4.2. 



We can see that Algorithm 2.1 could give a much shorter CP-decomposition 
than that given in (3.2). 

Example 4.3. Consider the nonnegativc P-matrix with diagonals partially given: 



(4.3) 



6 4 12 2 
4*013 
10 3 12 
2 11*1 
2 3 2 1 * 
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By Theorem 2.5, the corresponding is not A-full in this example. However, 

we still could obtain the following CP-completion matrix by using Algorithm 2.1 



6 


4 


1 


2 


2 




4 


4.1241 





1 


3 


5 


1 





3 


1 


2 




2 


1 


1 


1.0108 


1 


i=l 


2 


3 


2 


1 


5.2354 





with Ui and pi shown in Table 3. 



i 


Ui 


Pi 


1 


(0.4125, 0.0000, 0.3815, 0.2048, 0.0010) T 


6.2038 


2 


(0.0000, 0.0000, 0.4650, 0.1020, 0.4329) T 


9.6328 


3 


(0.4430, 0.3220, 0.0000, 0.1299, 0.1050) T 


23.9046 


4 


(0.1331, 0.3717, 0.0000, 0.0000, 0.4951) T 


11.9082 


5 


(0.1533, 0.0000, 0.0890, 0.3787, 0.3788) T 


1.7209 



Table 3. The points Ui and weights pi for A in Example 4.3. 



This also verifies the property of Algorithm 2.1 that if an E- matrix could be 
completed to a CP-matrix, no matter whether the corresponding M.[x]s is A-full or 
not, a CP-completion could always be obtained by Algorithm 2.1. 

Example 4.4. Consider the nonnegative P-matrix with one diagonal and all the 
off-diagonals given: 

" * 7 1 3 9 
7*585 

(4.4) 15*22 

3 8 2 3 1 
9 5 2 1 * 

By Theorem 3.3, (4.4) could be completed to a CP-matrix. We can also apply 
Algorithm 2.1 to obtain a CP-completion and a CP-decomposition as follows: 



12.0587 
7 
1 
3 
9 



7 

26.8175 
5 



1 

5 

4.3096 
2 
2 



9 
5 
2 
1 

11.5011 



5 

^2piU t uf, 
i=l 



where Ui and pi are listed in Table 4. 

Example 4.5. Consider the nonnegative P-matrix (cf. [1, Example 1.35]): 





(4.5) 



It was shown in [1] that (4.5) is not CP-completable. Since all the diagonals are 
given, we can use Algorithm 2.1 to verify this fact. Algorithm 2.1 terminates at 
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i 




Pi 


1 


(0.0017, 0.4932, 0.2021, 0.0988, 0.2042) T 


48.4520 


2 


(0.2505, 0.5393, 0.0062, 0.2040, 0.0000) T 


50.8482 


3 


(0.1796, 0.4880, 0.0000, 0.0917, 0.2407) T 


1.0233 


4 


(0.2154, 0.0000, 0.5541, 0.2305, 0.0000) T 


7.5805 


5 


(0.4869, 0.0000, 0.0000, 0.0000, 0.5131) T 


35.7825 



Table 4. The points and weights pi for A in Example 4.4. 



Step 2 with k = 3 as (2.10) is infeasible, which indicates that (4.5) is not CP- 
completable. 

Example 4.6. Consider the nonncgative P-matrix with the diagonals partially 
given: 



(4.6) 



The maximal principle submatrix of (4.6) is 



(4.7) 



1 1 2 

1 1 3 

2 3 3 

* * 3 



Since all the diagonals of (4.7) are given, we can use Algorithm 2.1 to verify whether 
(4.7) is CP-completable or not. Algorithm 2.1 terminates at Step 2 with k = 2, as 
(2.10) is infeasible. Thus (4.7) is not CP-completable, which indicates that (4.6) is 
not CP-completable cither by Proposition 3.3. 

Example 4.7. Consider the nonnegative E- matrix (3.3) which was mentioned in 
subsection 3.2. We rewrite it as follows: 



(4.8) 



1 1 2 

1 1 3 

2 3 * 



Though its maximal principal submatrix 

A' = 



is CP-completable, (4.8) is not. We considered the approximations of (4i 
form (3.5). Choosing Sh = 10~ k , we have the P-matrix sequence 



the 



(4.9) 



10~ fe 

1 

2 



1 

10~ fc 
3 



fc = l,2,. 
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By the proof of Theorem 3.6, we know (4.9) is CP-complctablc for all k. We used 
Algorithm 2.1 and obtained the following CP-completion matrices 

1.1 1 2 
Ax 1 1.1 3 

2 3 10.9524 



\iUiuf, 



A 2 = 



1.01 
1 

2 



1 2 
1.01 3 
3 56.2189 



and 



1=1 
2 

i=l 
4 

= ^2 <?iUiU?, 



T 



1.001 1 2 
A 3 = 1 1.001 3 

2 3 487.2967 
for fc = 1,2 and 3, respectively. The support points and the corresponding A- 
representing measure weights are shown in Table 5 and 6. 



i 


Ui 


A; 




Pi 


1 


(0.1254, 0.1881, 0.6866) T 


23.2350 


(0.0327, 0.0490, 0.9183) T 


66.6636 


2 


(0.6190, 0.3810, 0.0000) T 


1.9174 


(0.5124, 0.4876, 0.0000) T 


3.5753 



Table 5. The points and weights for A\ and A2 in Example 4.7. 



i 






1 


(0.5012, 0.4987, 0.0001) T 


3.7527 


2 


(0.0682, 0.0696, 0.8623) T 


11.6502 


3 


(0.0027, 0.0048, 0.9925) T 


485.9000 


4 


(0.4545, 0.3039, 0.2487) T 


0.0023 



Table 6. The support points and the weights for A3 in Example 4.7. 



Remark 4.8. For the case that the P-matrix have diagonals ungiven, if the maximal 
principle submatrix is not CP-complctable, then the P-matrix is not either. But 
even if the maximal principle submatrix is CP-completable, the P-matrix may not 
be CP-completable. Anyway, we know from Theorem 3.6 and Remark 3.7 that 
the P-matrix is in the closure of CP-completable P-matrices; there always exists a 
sequence of CP-complctablc P-matrices with the identifying vectors converging to 
that of the P-matrix. Moreover, Algorithm 2.1 could give not only a CP-completion 
but also a CP- factorization of the approximate CP-matrix. 

5. Conclusion 

In this paper, we have successfully solved the CP-completion problem for the 
nonnegative P-matrix. We formulate it as an £-TAMP and present an SDP Al- 
gorithm 2.1 for it. For the case that all the diagonal entries are given, since the 
corresponding M.[x]s is A-full, Algorithm 2.1 can give a certificate for the non- 
CP-completeness of the P-matrix. Moreover, if the P-matrix is CP-completable, 
Algorithm 2.1 is almost always guaranteed to get not only a CP-completion but 
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also a CP-decomposition. For the P-matrix with all diagonals ungiven, we show 
it can always be completed to a CP-matrix and Algorithm 2.1 can almost always 
obtain a CP-completion and a CP-decomposition with a smaller length. For the 
P-matrix with diagonals partially given, if it is not partial CP-completable, then a 
certificate can be obtained; if it is partial CP-complctable, then there always exists 
a sequence of CP-complctable P-matrices whose identifying vector converges to 
the identifying vector of the E- matrix, Algorithm 2.1 is almost always guaranteed 
to get a CP-completion if it is CP-completable. In particular, if the P-matrix has 
only one positive diagonal given, it is always CP-completable. 
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